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Abstract: Design of a doubly-clamped beam structure capable of 
localizing mechanical and optical energy at the nanoscale is presented. The 
optical design is based upon photonic crystal concepts in which patterning 
of a nanoscale-cross-section beam can result in strong optical localization 
to an effective optical mode volume of 0.2 cubic wavelengths ((A, c ) 3 ). By 
placing two identical nanobeams within the near field of each other, strong 
optomechanical coupling can be realized for differential motion between the 
beams. Current designs for thin film silicon nitride beams at a wavelength 
of X = 1.5 fjm indicate that such structures can simultaneously realize an 
optical g-factor of 7 x 10 6 , motional mass m u ~ 40 picograms, mechanical 
mode frequency Q.m/2k ~ 170 MHz, and an optomechanical coupling 
factor (goM = dC0 e /ch; = C0 c /Lom) with effective length Lom ~ A, = 1 .5 fim. 
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1. Introduction 

At a macroscopic level, the interaction of light with the mechanical degrees of freedom of a 
dielectric object can be calculated by considering the flux of momentum into or out of an object 
using the Maxwell stress-energy tensor. At a microscopic level, as in the case of atomic physics, 
one can define an interaction Hamiltonian between an atom and the light field in order to derive 
the various mechanical forces on the atom's center of mass, which in general depends upon 
both the external and internal degrees of freedom of the atom[lJ. In the case of a dielectric 
mechanical resonator, a direct relationship between the macroscopic dielectric and microscopic 
atomic theories can be made, and useful analogies may be forged|2] [3). The interactions of 
light with mechanically resonant objects is currently being actively explored in the field of 
cavity optomechanics as a means to obtain ground-state cooling of a macroscopic mechanical 
resonator l3l l4l l5l l6l 171 l8l l9l [lOl . The strength of optomechanical interactions in these system can 
be quantified on a per-photon basis by the rate of change of the cavity resonance frequency (co c ) 
with mechanical displacement amplitude (u), goM = d(O c /du = CO c /Lom- ^om is an effective 
length over which a cavity photon's momentum can be exchanged with the mechanical system. 
In this work we describe a simple doubly clamped nanobeam system (a so-called zipper cavity) 
which allows for the combined localization of optical and mechanical energy in a nanoscale 
structure so as to provide extremely large optomechanical coupling due to the gradient optical 
force. Optical energy is localized within the center of the cantilever using a one-dimensional 
photonic crystal in combination with total internal reflection. Beyond the analysis provided 
here, in the future, optimization of both the optical and mechanical properties of these chip- 
based structures should allow for a variety of new applications from precision metrology lfTTl to 
tunable photonics |[T2l IT3l . 

The outline for the paper is as follows. We begin with the optical design of a one-dimensional 
photonic crystal in a siliocn nitride nanobeam. Finite-element-method electromagnetic simula- 
tions are used to deduce the level of optical localization and the relevant optical losses within 
the struture. The mechanical properties of the zipper cavity are studied next, with numerical 
simulations used to determine the lower-lying mechanical eigenmodes. The tuning properties 
of a double nanobeam photonic crystal are then computed to estimate the strength of the op- 
tomechanical coupling for the differential in-plane motion of the beams. We conclude with a 
comparison of the zipper cavity properties with other more macroscopic optomechanical sys- 
tems, and a discussion of the future prospects for these sorts of chip-based gradient optical force 
devices. 

2. Optical design and simulation 

The optical design of the zipper cavity utilizes a quasi- ID photonic crystal structure to localize 
optical modes to the center of a nanoscale cross-section beam. There are a number of different, 
but related, design methodologies used to form low-loss optical resonances in such photonic 
crystal "nanowires" lfl4l [T5l [161 [171 [181 [T9l . In this paper we utilize concepts based upon an en- 
velope picture of the guided photonic crystal modes {20 1 in the patterned nanobeam. We choose 




Fig. 1. Bandstructure properties of the photonic crystal nanowire. (a) Axial bandstructure 
of the single beam photonic crystal structure, with nominal width and thickness, k is 
the wavevector of light in the direction of the ID photonic lattice, a is the lattice period, 
and ci/Xq is the "normalized" optical frequency. The light cone, denoted by the grey area 
and deliniated by the black light line, represents regions of frequency-wave-vector space 
in the bandstructure diagram in which light can radiate into the two transverse directions 
orthogonal to the axis of the photonic lattice. The two inset images show the electric field 
energy density of the valence (i) and conduction (ii) band-edge modes (the white outline is 
a contour plot of the refractive index of the nanowire). (b) Schematic of the double beam 
zipper cavity indicating the slot gap (s), the lattice constant (a), the beam width (w), and the 
axial and transverse hole lengths, h x and by, respectively, (c) Lattice constant (normalized 
to the lattice period in the mirror section of the cavity, a m ) versus hole number («;,) within 
the photonic lattice of the cavity. (c,d) Resulting frequency of the valence band-edge mode 
versus hole number. In (d) the normalized frequency in terms of the local lattice constant, 
(a/Xo) nii , is displayed. In (e) the local band-edge frequency is referenced to the valence 
band-edge in the mirror section of the cavity. The solid blue (dashed red) curve is the 
valence (conduction) band-edge. 



to work with "acceptor"-type modes formed from the lower-lying band of modes (the so-called 
"valence" band in reference to electronic semiconductor crystals) at the Brillouin zone bound- 
ary. In this way, the mode frequencies are as far away from the light line as possible, reducing 
leakage into the surrounding air-cladding of the nanobeam. Owing to the optically thin pho- 
tonic wires considered in this work, the modes of predominantly in-plane polarization (TE-like 
modes) are most strongly guided and are the modes of primary interest here. In order to form 
acceptor modes, the bandedge of the valence band must be increased near the center of the 
cavity, and decreased in the mirror sections defining the end of the cavity. This is a result of the 



negative group velocity dispersion of the valence band modes l20l . 

We have chosen to perform designs based upon thin films of silicon nitride, as opposed 
to higher refractive index materials such as silicon-on-insulator (SOI), for several reasons. One 
reason is that silicon nitride can be grown on silicon wafers with very high optical quality across 
a wide range of wavelengths covering the visible to the mid-infrared. We have measured lED 
optical (2-factors in excess of 3 x 10 6 for whispering-gallery modes of microdisks formed from 
stoichiometric silicon nitride deposited by low-pressure-chemical-vapor-deposition (LPCVD). 
A second reason is that LPCVD-deposited stoichiometric silicon nitride films on silicon have a 
large internal tensile stress, which has been shown to be critical in producing high-Q mechani- 
cal resonances in doubly-clamped nanobeams[22, 23 1. An additional concern is the two-photon 
absorption present in smaller bandgap semiconductors such as silicon and gallium arsenide, 
which results in additional free-carrier absorption (FCA), and which can result in significant 
parasitic effects such as thermo-optic and free-carrier dispersion in small-volume photonic crys- 
tal nanocavities[24|. Silicon nitride, with its large bangap (~ 3 eV), requires three (as opposed 
to two) 1 eV photons to be absorbed simultaneously, greatly reducing nonlinear absorption in 
the near-IR. An obvious drawback of using silicon nitride thin films, is the lower refractive in- 
dex (n ~ 2) of these films in comparison to semiconductor films (« ~ 3.4). As is shown below 
(and in Ref. |[T9l ). with carefully chosen designs, high-g photonic crystal optical cavities of 
sub-cubic-wavelength mode volume can still be formed in silicon nitride thin films. 

The bandstructure of a single beam silicon nitride nanowire, calculated using the MIT pho- 
tonic bands software package|25], is displayed in Fig.Q] As described in more detail below, 
the simulation is performed for a "nominal" structure defined by a lattice normalized beam 
thickness (t = t/a = 2/3) and beam width (VP = w/a = 7/6). The refractive index of the silicon 
nitride beam is taken as n = 2, and a resolution of 64 points per axial lattice period is used to 
ensure accurate band frequencies. Our coordinate convention is (see Fig. [3}: 00 x the in-plane 
coordinate along the long axis of the cavity, (ii) y the in-plane transverse coordinate, and (iii) z 
the out-of-plane transverse coordinate. Only the lower-lying bands with modes of even parity 
in the ^-direction and odd parity in the y-direction are shown, corresponding to the fundamental 
TE-like modes of the beam waveguide. As indicated by the electric field energy density plots of 
the two lowest lying band-edge modes (inset to Fig. da)), the valence band lies predominantly 
in the region of the high refractive index silicon nitride beam, whereas the upper "conduction" 
band mode lies predominantly in the region of the air hole patterning. 

The formation of localized optical cavity resonances is accomplished by introducing a "de- 
fect" into the photonic lattice. The defect region in the structures studied here consists of a 
quadratic grade in the lattice constant of the linear array of air holes near the cavity center. 
In order to reduce transverse radiation loss, we choose to use a defect which supports an odd 
symmetry fundamental mode along the axial direction. Since the valence band modes tend to 
have electric field intensity predominantly inside the high-dielectric region (and nodes of the 
electric field in the low-dielectric constant air holes), a defect in which an air hole is at the 
center of the cavity yields an odd parity fundamental mode along the axial direction. Here, and 
in what follows, we use a cavity defect region consisting of the central 15 holes, with the lattice 
period varied from a nominal value in the outer mirror section (a m ) to 90% of the nominal 
value at the center of the defect region. This was found to provide a good balance between 
axial-localization of the cavity modes and radiation loss into the y-z transverse directions of 
the nanobeams. In Fig. 02c) we plot the local lattice period, defined as a„ h = x(ni, + 1) — x{tih), 
versus air hole number tit, along the length of the cavity. In Fig. Old) we plot the corresponding 
lattice-normalized frequency of the local TE-like valence band-edge modes of the single beam 
cavity. The small variation in lattice-normalized frequency is a result of the distortion in the 
aspect ratio of the (perfectly periodic) structure as the lattice period is changed. More useful 



is the plot in Fig. [TJe) which shows the frequency of the local valence band-edge (blue solid 
curve) and conduction band-edge (red dashed curve) modes normalized to the valence band- 
edge mode frequency in the mirror section of the cavity. The quadratic grade in lattice constant 
results in a nearly-harmonic shift in the valence band-edge frequency versus position in the cen- 
ter of the cavity, with the band-edge frequency at the cavity center lying approximately mid-gap 
between the valence and conduction band-edges in the mirror section of the cavity. 




Fig. 2. Optical design principle of the zipper cavity, (a) Axial bandstructure of the double 
beam quasi- ID photonic crystal structure, with nominal width, thickness, and slot gap. The 
blue curves are the bonded bands and the red curves are the anti-bonded bands. Valence and 
conduction band-edge modes of the bonded (b,c) and anti-bonded (d,e) bands, respectively, 
(f) Illustration of the defect cavity formation at the Brillouin-zone boundary. The splitting 
between the two manifolds is indicated by AX + ,_. 

For the double-beam design of the zipper cavity, two photonic crystal nanowires are placed 
in the near-field of each other as illustrated in Fig. |TJb). The strong coupling between two 
nearly-identical nanobeams results in a bandstructure consisting of even and odd parity super- 
positions of the TE-polarized single beam photonic bands (Fig. |2ja)). We term the even parity 
supermodes, bonded modes, and the odd parity supermodes, anti-bonded modes fPH . The anti- 
bonded manifold of resonant modes are shifted to higher frequency than the bonded manifold 
of modes, with splitting (AA, +j _) being dependent upon the slot gap (s) between the nanobeams. 
The electric field energy density plots of the bonded valence and conduction band-edge modes 
are shown in Figs. |2jb) and |2jc), respectively. The corresponding anti-bonded modes, with 
noticeably reduced energy density in the slot gap, are displayed in Fig. |2td,e). Owing to the 
negative curvature of the valence band forming the localized cavity modes, the fundamental 
mode of the cavity for each manifold has the highest frequency, with the higher-order modes 
having reduced frequencies (schematically illustrated in Fig. |2jf)). In what follows, we will 
primarily be interested in the bonded mode manifold due to the larger electric field intensity of 
these modes in the slot gap. 

With the zipper cavity now defined, we perform a series of finite-element-method (FEM), 
fully vectorial, 3D simulations of the localized cavity modes ll26l . In these simulations we use 
the graded cavity design described above. The numerical mesh density is adjusted to obtain 
convergent values for the frequency and optical Q-factor of the cavity modes, and the structure 
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Fig. 3. Transverse electric field (E v ) mode profile of the fundamental bonded mode (TE + q): 
(a) top- view, (b) cross-section. Transverse electric field (E v ) mode profile of the fundamen- 
tal anti-bonded mode (TE_,o) : ( c ) top-view, (d) cross-section. The field colormap corre- 
sponds to +1 (red), (white), and -1 (blue). Electric field energy density of the (e) fun- 
damental (TE + o), (f) second-order (TE + i), and (g) third-order (TE + 2) bonded optical 
modes. The intensity colormap ranges from +1 (red) to (blue). 



is simulated with all available symmetries taken into account allowing for l/8th the simulation 
volume. Scattering boundary conditions are used in the outer axial and transverse boundaries 
to provide a nearly-reflectionless boundary for out-going radiation. In the transverse direction 
(y-z) a cylindrical outer boundary is also used to reduce the total simulation volume. Finally, a 
check of the accuracy of our FEM simulations was also performed through a series of equivalent 
simulations using finite-difference time-domain code (Lumerical |2"7| ). Figures [3ja) and 0c) 
display the FEM-calculated electric field mode profiles of the fundamental bonded (TE + 0) 
and anti-bonded (TE_,o) modes, respectively, of the double beam zipper cavity. Cross-sectional 
electric field profiles, displayed in Figs. [3jb,d), clearly show the even and odd parity of the 
modes. Also shown in Figs. Oe-g) are the electric field intensity of the lowest three bonded 
cavity mode orders, TE+,0, TE+ i i, and TE+,2. 

Design variations of the zipper cavity are performed around a "nominal" structure with the 
following (normalized) dimenensions: (i) beam width, w = w/a m = 700/600, (ii) beam thick- 
ness, t = t/a m = 400/600, (iii) axial hole length, h x = h x /a m = 267/600, (iv) axial hole width, 
h y = h y /a m = 400/600, and (v) slot gap, s = s/a m = 100/600, where a m is the lattice periodic- 
ity in the cavity mirror section. The length of the beam is set by the number of air hole periods 
in the cavity, which for the nominal structure is N/, = 47 (23 holes to the left and right of the 
central hole, with the central 15 holes defining the defect region). The filling fraction of the air 
holes in the nominal structure is / = 25.4%. For a fundmamental bonded mode wavelength of 



X ~ 1500 nm, the nominal lattice constant is a„, = 600 nm, hence the normalization by 600 in 
the above expressions for normalized dimensions. 

Figure H shows the simulated £>-factor of the fundamental bonded cavity mode (TE+,o) ver- 
sus the total number of periods of the cavity for the nominal structure. The normalized 
cavity resonance frequency is calculated to be a m jX c = 0.3927. The radiation loss from the 
cavity is broken into two parts, yielding two effective Q-factors: (i) the axial radiation loss out 
the ends of the nanobeams (yielding Qn), and (ii) the radiation loss transverse to the long axis 
of the zipper cavity and intercepted by the transverse boundary (yielding Q±). In Figure|H £>| 
is seen to rise exponentially as a function of Nh, with an order of magnitude increase in g-factor 
for every 6 additional periods of air holes. The transverse Q is seen to rise initially with hole 
periods, but then levels off and saturates at a value of Qj_ = 7x 10 6 . The variation in Qj_ for 
structures with small Nh, and low Qn, is a result of weak coupling between radiation loss into 
these (arbitrarily chosen) directions. Small reflections at the end of the nanobeams results in 
a small amount of axial radiation making its way to the transverse boundary. Nevertheless, a 
structure with N/, > 47 results in a total radiation Q-factor limited by the transverse Q-factor of 
7 x 10 6 . This value is very large given the modest refractive index of the silicon nitride film and 
small cross-section of the nanobeams (large air filling fraction). 
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Fig. 4. Optical g-factor (axial, transverse, and total) versus number of hole periods in the 
cavity, N/,. The nominal structure corresponds to the maximum hole number in this plot, 
N h = 47. 




In order to study the dependence of g-factor on the hole size, we have also simulated the 
nominal structure with varying axial and transverse hole size, as shown in Figs. [5] and [6] In 
each of the plots the axial Q-f&clor: increases with increasing hole size, but then saturates, as 
the hole size approaches that of the nominal structure. For hole sizes larger than the nominal 
structure, the g| slightly drops, as does the transverse Q-factor. This drop in transverse g-factor 
is a result of the increased normalized frequency of the mode (higher air filling fraction), which 
pushes the mode closer to the air-cladding light line, increasing radiation into the cladding 
(the drop in Q\\ is a result of low-angle transverse radiation making its way to the boundary at 
the end of the nanobeams). Therefore, the nominal structure represents an optimal structure in 
this regard, allowing for tight axial localization of the mode, without decreasing the transverse 
g-factor. 

The strength of light-matter interactions depend upon a position-dependent effective optical 
mode volume of a resonant cavity, 

J £ (r)|E(r)| 2 d 3 r 

Veff(l"0) = -, 77=-, TTj , (1) 

e(r )|E(r )| 2 
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Fig. 6. (a) Normalized frequency and (b) Optical g-factor (axial, transverse, and total) 
versus normalized transverse hole length, h y . The nominal structure is indicated by the 
dashed green line. 



where e(r) is the position dependent dielectric constant and ro is the position of interest. In the 
case of the zipper cavity, the strength of the optomechanical coupling depends upon V e ff(ro) 
evaluated in the slot gap of the nanobeams. Similarly, in the field of cavity-QED, the effective 
mode volume can be used to estimate the coherent coupling rate between an "atom" and the 
cavity field. In Fig. [7] we plot the effective mode volume versus slot gap size, s, between the 
nanobeams for the TEo,+ mode of the nominal structure. We plot two different effective mode 
volumes: (i) V g , the effective mode volume evaluated at the center of the nanobeam gap near 
the center of the cavity where the field is most intense, and (ii) V p , the minimum effective 
mode volume evaluated at the position of peak electric field energy density in the cavity. For 
normalized gap widths, s < 1/6 the modes of the nanobeam are strongly coupled resulting in 
a peak electric field intensity in between the slot gap (hence the two effective mode volumes 
track each other). The effective mode volumes approach a value of V e g = 0.1 (A, e ) 3 as the slot gap 
approaches zero. This small value is a result of the discontinuity in the dominant polarization 
of the TE mode (E y ) as it crosses into the slot gap l28l . For slot gaps s > 1/6, the minimum 



effective mode volume saturates at a value of V p = 0.225 (A. c ) 3 corresponding to that of a single 
nanobeam (i.e., no enhancement from energy being pushed into the slot gap). The effective 
mode volume evaluated at the center of the slot gap, on the otherhand, continues to rise with 
slot gap due to the exponential decay of the field in the gap. 
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Fig. 7. Effective mode volume versus normalized slot gap, s. The nominal structure is 
indicated by the dashed green line. 

A final variation considered is the nanobeam width. In Fig.|8]the nominal structure is varied 
by adjusting the normalized beam width while holding the filling fraction of the air holes fixed. 
This is done by scaling the the transverse hole length with the beam width. The axial g-factor 
increases with the beam width due to the increased effective index of the guided mode, and 
thus increased effecive contrast of the quasi-lD photonic crystal. Due to the reduced lateral 
localization of the cavity mode, the effective mode volumes also increase with beam width. In 
contrast, the transverse g-factor remains approximately constant. 
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Fig. 8. (a) Normalized frequency, (b) g-factor and (c) Effective mode volume versus nor- 
malized beam width. The nominal structure is indicated by the dashed green line. 



3. Mechanical mode analysis 



The mechanical modes of the zipper cavity can be categorized into common and differential 
modes of in-plane (labeled h) and out-of-plane motion (labeled v) of the two nanobeams. In 
addition there are compression (labeled c) and twisting (labeled t) modes of the beams. In this 
work we focus on the in-plane differential modes, h q d, as these modes are the most strongly 
opto-mechanically coupled due to the large change in the slot gap per unit (strain) energy. 
We use FEM numerical simulations to calculate the mechanical mode patterns and mode fre- 
quencies, the first few orders of which are shown in Fig. [9] The material properties of silicon 
nitride for the FEM simulations were obtained from a number of references. Where possible 
we have used parameters most closely associate with stoichiometric, low-pressure-chemical- 
vapor-deposition (LPCVD) silicon nitride deposited on (100) Si: mass density p = 3100 kg/m 3 , 
Young's modulus Y ~ 290 GPa, internal tensile stress a ~ 1 GPa[22 23 1, coefficient of ther- 
mal expansion t\te =3.3 x 10~ 6 K \ thermal conductivity K,/, ~ 20 W/m/K, and specific heat 



For the mechanical mode properties tabulated in Table Q] we have analyzed the nominal zip- 
per cavity structure at an operating wavelength of X 1.5 /Jin, corresponding to a geometry 
with a m = 600 nm, t = 400 nm, s = 100 nm, w = 700 nm, h x = 267 nm, and h y = 400 nm. 
The total number of air holes is set to Nj, = 55, ensuring a theoretical Q-factor dominated by 
transverse radiation (Q±), yielding a cavity length of / = 36 /jm. The zipper cavity is clamped 
at both ends using a fixed boundary conditions at the far ends of the "clamping pads" seen in 
Fig- El This clamping scheme is suitable for estimating the mechanical mode eigenfrequencies, 
although more complex clamping schemes envisioned for real devices will likely introduce 
modified splittings betweeen nearly-degenrate common and differential modes. The effective 
spring constant listed in Table Q] for each mode is based upon a motional mass equal to that 
of the true physical mass of the patterned nanobeams, m u « 43 picograms (see below for self- 
consistent definition of motional mass). 

The resulting frequency for the fundamental h\ c i mode is £1/2% w 8 MHz. The in-plane mode 
frequency of a doubly clamped beam, with I 3> w,t, is approximately given by ||22l : 



where q is the mode index (approximately an integer), A = tw is the cross-sectional area of the 
beams, and L = fw 3 /12 is the cross-sectional moment of inertial about the in-plane axis (y) of 
the beam. This fits the numerical data reasonably accurately assuming an effective beam width 
of w' = (1 —f)w. From the scaling in eq. (|2), one finds that for mode number q > 3 (where 

internal stress can be neglected and Q, q /2n w i/T^/js") that the in-plane frequency scales in- 
versely with the square of the beam length, linearly with the beam width, and independent of 
the beam thickness. Therefore, a linear increase in the resonant frequency can be obtained by 
moving to shorter optical wavelengths and scaling the structure with wavelength. For much 
larger increases in mechanical resonance frequency, one must resort to higher-order in-plane 
modes of motion. Optomechanical coupling to these modes is discussed below. 

We have also studied the expected thermal properties of the zipper cavity, again assuming a 
1 .5 /jm wavelength of operation. Due to the air-filling-fraction of the etched holes in the zipper 
cavity nanobeams, the thermal conductivity of the patterned beams is approximately F t i, = 75% 
of the bulk value. A simple estimate for the thermal resistance of the zipper cavity is R t h ~ 
I / (StwTthKth) ~ 1-15 x 10 6 K/W, where the factor of 1/8 comes from the ability for heat to 
escape out either end of the nanobeams and in either direction. The physical mass of the zipper 
cavity, taking into account the etched holes, is approximately m = 43 picograms. The heat 



c sh = 0.7 J/g/K. 




(2) 




Fig. 9. Mechanical eigenmode displacement plots: (a) lst-order in-plane common mode, 
(b) lst-order in-plane differential mode, (c) lst-order out-of-plane common mode, and (d) 
lst-order out-of-plane differential mode. The color represents total displacement ampli- 
tude, as does the deformation of the structure. The arrows indicate the local direction of 
displacement. 



capacity of the zipper cavity is then q,«3x 10 J/K. From the heat capacity and the thermal 
resistance, the thermal decay rate is estimated to be y,/, = l/R t f,Cf, ~ 2.9 x 10 4 s _1 . Finite- 
element-method simulations of the thermal properties of the zipper cavity yield an effective 
thermal resistance of R t h = 1.09 x 10 6 KAV and a thermal decay rate of y,/, = 5.26 x 10 4 s _1 
for temperature at the center of the zipper cavity, in reasonable correspondence to the estimated 
values. 

4. Optomechanical coupling 

With the optical and mechanical modes of the zipper cavity now characterized, we proceed to 
consider the optomechanical coupling of the optical and mechanical degress of freedom. As 
described at the outset, the parameter describing the strength of optomechanical coupling is the 
frequency shift in the cavity mode frequency versus mechanical displacement, goivt = d(O c /dw, 
where <B t - is the cavity resonance frequency and u represents an amplitude of the mechanical 
displacement. In the case of the commonly studied Fabry-Perot cavity 181. goM = 03 C /L C , where 
L c is approximately the physical length of the cavity. A similar relation holds for whispering 
gallery structures, such as the recently studied microtoroid||2), in which the optomechanical 
coupling is proportional to the inverse of the radius of cavity (R), goM = d> c /R. Both these 
devices utilize the radiation pressurce, or scattering, force of light. By comparison, the zipper 
cavity operates using the gradient force for which the optomechanical coupling length can be 
on the scale of the wavelength of light, Lom ~ ^e- Similar to the scaling found in the previous 
section for the effective mode volume, the optomechanical coupling length scales exponentially 
with the slot gap, Lom ~ ^-c<? as , where a is proportional to the refractive index contrast between 
the nanobeams forming the zipper cavity and the surrounding cladding. 

In Fig. [Tol a) we plot the tuning curve for the nominal zipper cavity structure studied in the 
previous two sections versus the normalized slot gap width, s = s/ a m . Due to the strong intensity 
of the bonded mode in the center of the slot gap, it tunes more quickly than the anti-bonded 
mode. This tuning curve can be used to estimate the optomechanical coupling for the in-plane 



Table 1. Summary of mechanical mode properties. 



Mode label 


Q. M /2n (MHz) 


k u (N/m) 


O *_J1V1 \ / * W C / 


1 1 i 


5.95 


60 




Vij 


6.15 


64 




1 _t 


12.3 


257 




VTri 


12.7 


274 




1 Jl 


19.2 


626 




V3d 


20.0 


679 




"-Yd 


7.91 


106 


1.24 


h\ c 


7.94 


107 


-0 


hid 


18.2 


562 


-0 


fl2c 


18.3 


568 


-0 


hd 


31.8 


1717 


1.16 


he 


32.0 


1738 


-0 


h9d 


167.7 


4.77 x 10 4 


0.63 


hg c 


168.0 


4.79 x 10 4 


-0 


hd 


41.0 


2854 




he 


41.1 


2868 




Clc 


78.6 


1.05 x 10 4 




C\d 


79.4 


1.07 x 10 4 





differential mode of motion of the zipper cavity nanobeams. The in-plane common mechanical 
modes and both types of vertical mechanical modes are expected to provide a much smaller 
level of optomechanical coupling due to the reduced change in slot gap with these types of 
motion. For complex geometries and motional patterns, one must use a consistent definition of 
displacement amplitude, u, in determing goM, '«eff (motional mass), and A: e ff (effective spring 
constant). In this work we use a convention in which u(t) represents the amplitude of motion 
for a normalized mechanical eigenmode displacement field pattern: 

f (x) 

u„{r,t) =u n (t) (3) 
V}/ '|f„WI 2 d* 

where n is a mode label, / is the length of nanobeam, and, for the simple beam geometry 
considered here, the displacement vector is only a function of position along the long axis of 
the nanobeams (x). With this definition of amplitude, the effective motional mass is simply 
the total mass of the two nanobeams (m u = 43 picograms), and the effective spring constant is 
defined by the usual relation k u = m u Dr u , with Q.m the mechanical eigenmode frequency. The 
amplitude associated with zero-point motion and used in the equipartition theorem to determine 
the thermal excitation of the mechanical mode is then u n (t). Note, an alternative, but equally 
effective method, defines the amplitude first, and then adjusts the effective motional mass based 
upon the strain energy of the mechanical motion. 

Our chosen normalization prescription yields (approximately) for the gth odd-order in- 
plane differential mechanical mode, u° 2 d (z,t) « Uh qi {t) (y\ cos(qttx/l) + y2Cos(qnx/l)), where 
y\ and yi are (transverse) in-plane unit vectors associated with first and second nanobeams, 



respectively, and which point in opposite directions away from the center of the gap between 
the nanobeams. The even-order modes are anti-symmetric about the long axis of the cavity 
and are given approximately by, u| d (z,t) « u h qi {f) ($1 sin(qKx/l) + y2&in(qnx/l)). To be con- 
sistent then, with this definition of mode displacement amplitude, goM must be defined in 
terms of the rate of change of cavity frequency with respect to half the change in slot gap 
(gOM ~ dcOc/^dSs), as the amplitude u% d (t) corresponds to a (peak) change in slot gap of 

Fig. fTOl b) plots the optomechanical coupling length for each of the bonded and anti-bonded 
fundamental modes from the derivative of their tuning curves in Fig. [Tol a). This plot shows 
that for a normalized slot gap of s = 0.1 (or roughly s = 0.04A, e ), the optomechanical coupling 
length to the fundamental bonded optical mode can be as small as Lo^/a m « 2. For the nor- 
malized frequency of the bonded mode (a m /X c w 0.4), this corresponds to Lom ~ 0.8^ c , as 
expected from the arguments laid out in the introduction. The TE_ o has a significantly smaller 
optomechanical coupling due to its reduced electric field energy in the slot. 

An estimate of the optomechanical coupling to the different in-plane differential mechani- 
cal modes can be approximated by averaging the displacement amplitude field pattern of the 
mechanical mode against the (normalized) optical intensity pattern of the zipper cavity op- 
tical modes|29|. Given the odd symmetry of the even-order li qc / mechanical modes, and the 
even symmetry of the optical intensity for the zipper cavity optical modes, the optomechanical 
coupling to the even-order h qc \ modes is approximately zero. The odd-order in-plane differ- 
ential modes, on the otherhand, have an anti-node of mechanical displacement at the optical 
cavity center and an even long-axis symmetry. For mode numbers small enough that the half- 
wavelength of the mechanical mode is roughly as large, or larger, than the effective length of 
the optical cavity mode along the axis of the beam, the optical mode will only sense the central 
half-wave displacement of the mechanical mode and the optomechanical coupling should still 
be quite large. As an example, from the intensity plot of the TE + o fundamental bonded optical 
mode in Fig. [3] the effective length of the optical mode along the long-axis of the nanobeams is 
L e s ~ 1a m = 4.2 /jm (for X c = 1.5 /jm). The mechanical mode index q is roughly equal to the 
number of half-wavelengths of the mechanical mode along the axis of the zipper cavity. There- 
fore, for the zipper cavity of length / = 36 /jm studied above, the 9th order in-plane differential 
mechanical mode with £lf u , i /2n w 170 MHz has a half-wavelength equal to 4 fim « L e ff- The 
resulting optomechanical coupling of the hgj mechanical mode to the TE + ,o optical mode is 
then still relatively large, equal to approximately half that of the coupling to the fundamental 
h\d mechanical mode. The optomechanical coupling factor to the TE+ ) o for each of the in-plane 
differential modes is tabulated along with the mechanical mode properties in Table Q] 

There are several physical ways of understanding the strength of the optomechanical cou- 
pling represented by goM- The per-photon mechanical force is given by F p h = /zgoM- For the 
zipper cavity, this yields a per-photon force of F p h f« h(O c /X c , which at near-infrared wave- 
lengths corresponds to 0.2 pN/photon. Such a force could be measured using other, non-optical 
techniques, and could provide a means for detecting single photons in a non-demolition man- 
ner. Also, through the optomechanical coupling, intra-cavity light can stiffen ll30l l7l [3Tl[32l and 
dampen|4l[5l|6]0[8]|9] the motion of the coupled mechanical oscillatorf2). A perturbative anal- 
ysis shows that in the sideband unresolved limit (Q.m <C T) the effective mechanical frequency 
(Cl' M ) and damping rate (7^) are given by the following relations (see Ref. J2j): 



Ym = Ym - 



2 l fl 0| gQM 

A 2 co c .m„ 

4|«o| 2 g^M r 
A 4 <o c m u 



A' , 



(4) 
(5) 



where £1« and Jm are the bare mechanical properties of the zipper cavity, \ao\ is the time- 
averaged stored optical cavity energy, A', = (0/ — co ( is the laser-cavity detuning, F is the 
waveguide-loaded optical cavity energy decay rate, and A 2 = (A' fJ ) 2 + (T/2) 2 . The maxmi- 
mum "optical spring" effect occurs at a detuning point of A' a = T/2. For this laser-cavity 
detuning, a single cavity photon introduces a shift in the mechanical frequency correspond- 
ing to A(£l 2 M )/Sl 2 M = (2Qho) c ) I \X 2 c k eS ). For the zipper cavity, with Q = 5 x 10 6 , X c = 1 /jm, 
and k/ lhl = 100 N/m, the resulting single-photon mechanical frequency shift is approximately 
A(Q.m)/Qm =1%- Thus, even a single cavity photon may yield a measurable shift in the me- 
chanical resonance frequency. 
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Fig. 10. (a) Bonded and anti-bonded mode tuning curves versus normalized nanobeam gap. 
(b) Corresponding normalized effective optomehcanical coupling length, Lom = i-OM/ fl m- 
The nominal structure is indicated by the dashed green line. 



5. Summary and Discussion 

Using photonic crystal concepts, we have designed an optomechanical system in which optical 
and mechanical energy can be co-localized in a cubic-micron volume and efficiently coupled 
through the gradient optical force. In the particular design studied here, a "zipper" cavity con- 
sisting of two nanoscale beams of silicon nitride, doubly clamped and patterned with a linear 
array of air holes, is used to form the optical cavity and the mechanical resonator. Mechanical 
oscillations consisting of differential motion of the doubly-clamped silicon nitride nanobeams 
results in an optomechanical coupling constant as large as goM ~ CD C /A, C -, where co c and X c are 
the optical resonant cavity frequency and wavelength, respectively. This coupling is several or- 
ders of magnitude larger than has been demonstrated in high-Finesse Fabry-Perot cavities, and 
is more than order of magnitude larger than for whispering-gallery micrototoid structures, both 
of which rely upon the radiation pressure force. Finite-element-method (FEM) simulations of 
the zipper cavity show that a structure with an optical Q = 5 x 10 6 , mechanical resonance fre- 
quency of Q.m/2k 170 MHz, and motional mass of m u ~ 40 picograms is possible. In the 
future, further increase in the mechanical frequency and reduction in the motional mass may be 



attained by using planar phononic crvstals ll33l to form the mechanical resonator. The combi- 
nation of phononic and photonic crystals would also provide an integrated, chip-scale platform 
for routing and coupling optical and mechanical energy. 

Beyond cavity optomechanics, the zipper cavity may also find application in the field of 
cavity QED. In particular, the zipper cavity as described here is suitable for a broad range of 
wavelengths from the visible to the mid-infrared. The optical mode volume is made smaller by 
the sub-wavelength slot gap between the nanobeams l28l . with V e ff ~ 0.2(X^) for a slot gap of 
s r-j X C /1Q. As an example, one can imagine placing nanoparticles of diamond (a popular solid- 
state system for quantum information processing |34l [35l [36l l37l [38l [39l ) in the gap between 
the nanobeams. Such "pick-and-place" techniques have been used with other, larger, optical 
cavities with good success l40ll4Tll42l . In the zipper cavity case, the small V c ff would produce a 
coherent coupling rate with the zero-phonon-line (ZPL) of the NV~ transition of approximately, 
^zpl/2ji ~ 3 GHz, even after accounting for the 3 — 5% branching ratio for the ZPL line. This is 
more than 100 times the radiative-limited linewidth measured for the NV~ transition (12 MHz), 
and more than 10 times the theoretical zipper cavity decay rate (90 MHz), putting the coupled 
system deep within the strong coupling regime. The additional benefit provided by the zipper 
cavity is the ability to rapidly tune the cavity frequency into and out of resonance with the ZPL 
of the NV~ transition. If mechanical resonance frequencies could be increased towards GHz 
values, using the suggested phononic crystal concepts for instance, then new approaches to 
photon-mediated quantum interactions and quantum state transfer can be envisioned for solid- 
state cavity QED systems. 
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